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ABSTRACT 

We report on the measurement of the 7.5-14.7 micron spectrum for the tran- 
siting extrasolar giant planet HD 189733b using the Infrared Spectrograph on 
the Spitzer Space Telescope. Though the observations comprise only 12 hours 
of telescope time, the continuum is well measured and has a flux ranging from 
0.6 mJy to 1.8 mJy over the wavelength range, or 0.49 ± 0.02% of the flux of 
the parent star. The variation in the measured fractional flux is very nearly flat 
over the entire wavelength range and shows no indication of significant absorp- 
tion by water or methane, in contrast with the predictions of most atmospheric 
models. Models with strong day/night differences appear to be disfavored by the 
data, suggesting that heat redistribution to the night side of the planet is highly 
efficient. 

Subject headings: Stars: — Planetary Systems, Stars: Binaries: Eclipsing, Stars: 
Individual (HD 189733) 



1. Introduction 



HD 189733 is a nearby KOV star with a magnitude of K = 5.5 and an expected flux 
of 250 mJy at 10 mic rons. HD 189733b w as initially discovered in the ELODIE search for 
transi ting exoplanets ( Bouchv et al. 20051) an d subsequently detected in existing Hipparcos 
data iHebrard fc Lecavelier Pes Etangs Ikood ). In addition to HP 189733b, HP 189733 is 
orbited by an M dwarf at a distance of 216 AU (IBakos et al. Il2006al ). HP 1897 33b orbits its 
parent star in 2.2 days and has an estimated radius of R p = 1.154 ± 0.033i?j (IBakos et al 
2006bl ). From first to last contact, the duration of the secondary eclipse is ~ 1.9 hours. 



The planet-to-star surface area ratio for HP 189733b is ~ 0.024 (IBakos et al. Il2006bl ). 



Owing to an orbital radius of 0.031 AU, the upper atmospheric temperature of the planet 
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is > 1000 K (IFortney et al. Il2006al ). Together these factors yield an expected planet-to-star 
flux ratio that is almost a factor of two greater at 10 microns than for either TrES-1 or HD 
209458b, two other exop lanets whose thermal emi s sion has recently be en detected by the 



Spitzer Space Telescop e (ICharbonneau et al. 1120051 ; iDeming et al. 1 120051 ) . This expectation 



has been borne out by Deming et al. ( 20061 ). who used the Peak- up Imaging array portion 
of Spitzer 's Infrared Spectrograph (IRS) to measure a secondary eclipse depth at 16/iin of 
0.551 ±0.030%. This makes HD 189733b a prime candidate for IRS spectroscopy and a more 
detailed characterization of its physical properties. 

We briefly describe the observations in Section [2] and the analysis in Section [3j We 
discuss the spectrum of HD 189733b in Section 0] 



2. Observations 

HD 189733 was observed with the Infrared Spectrograph (IRS) on October 21st and 
November 21st, 2006 as part of General Observer program #30473. In each case the star 
was observed for a six- hour span centered on the predicted time of secondary eclipse. All 
observations were made in first order (A = 7.4 — 14.5/Um) using the Short-Low module, giving 
a spectral resolution of ~ 80. The performance and capabilities of the Short-Low module, 
which has th e clean est and most sensitive of the IRS detectors, is discussed in detail by 



Houck et al. I fl2004h . 



The Oct 21st observations were carried out in "Staring" mode, wherein the telescope 
was periodically repointed to place the target alternately at two dither positions along the 
slit. Alternating between two slit positions provided some insurance against the possible 
large buildup of latent charge (see below). The observations were preceded by a high- 
accuracy peak-up approximately 3 hours before secondary eclipse and were taken continu- 
ously through the out-of-eclipse, ingress, secondary eclipse, egress, and out-of-eclipse portions 
of HD 189733b's orbit. Thirty successive 14-second exposures were taken at each slit position 
before moving to the next dither position. A total of 30 telescope moves were commanded, 
yielding 450 independent integrations at each dither position over the course of six hours. 

A preliminary analysis of the Oct 21st observations showed that charge latents were 
considerably more benign than initial estimates. An observing sequence scheduled for Nov 
21st was consequently modified to place the target at the center of the Short-Low slit for the 
entire observing sequence. Small pointing errors associated with any commanded telescope 
movement (which contribute significantly to the noise budget) were thus avoided. Having 
thereby also eliminated telescope slew and settle times, we were able to fit a total 950 14- 
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second exposures into the six hours available. 



3. Analysis 

Our analysis makes use of the two-dimensional Basic Calibrated Data products from the 
S14 version of the IRS pipeline. Using Version 1.4 of the SPICE spectral extraction tool, we 
apply the optimal extraction feature to extract one-dimensional spectra for each of the 1850 
exposures separately. We fix the center of each extraction window to take advantage of the 
pointing performance of the telescope and avoid adding noise associated with a measurement 
of the sub-pixel position of the star for each exposure. 

For reasons that will become clear, we treat the integrations at the two dither positions 
on Oct 21st and the single-pointing integrations taken on Nov 21st as three separate sets 
of observations. In principal, it should be possible to simply coadd all integrations in each 
observation set taken during eclipse (starlight only) and subtract them from a coaddition of 
all spectra taken out of eclipse (star plus planet). In practice, due to the limited integration 
time available when the planet is in eclipse, several systematics can manifest themselves as 
either increased noise or spectral offsets. 

The Spitzer Space Telescope's stationary pointing performance is not perfect. Due 
to a small, cyclic heat source within the spacecraft structure, there is a periodic varia- 
tion in alignment between the Star Tracker assembly and the telescope boresight. This 
l eads to a pointing oscillation w ith a period of about an hour and a magnitude of ~ 0.05" 



([Morales- Calderon et al. 1 120061 ) . Though small, this oscillation has significant consequences 
for observations that require precision on the order of 0.1%. In addition to variations in the 
total amount of light entering the slit, we are subject to significant signal modulation due 
to drifts across pixels with different (and poorly calibrated) sensitivities. 

We analyze each wavelength bin as an independent time series and then combine the 
results to produce a final ratio spectrum. As an example, we show in Figure 1 the measured 
flux as a function of time at dither position 1 for the wavelength bin centered on 7.93/im 
. Several characteristics that affect all wavelength bins to a greater or lesser degree are in 
evidence: 1) There is a rapid rise in the counts during the first 10 minutes of the sequence 
that is clearly related to latent charge buildup. 2) There is a longer term variation that is 
probably related to both the buildup and decay of latent charge, as well as to possible drifts 
in and out of the slit. 3) There is a cyclic behavior due to the pointing oscillations. 

We have found that the buildup and decay of latent charge can vary greatly from one 
pixel to the next. Rather than trying to model the ensemble, short-term latent behavior for 
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each wavelength bin, we simply exclude the first 10 minutes of data from our analysis. Since 
our final signal-to-noise ratio is determined primarily by data taken during secondary eclipse, 
the loss of the first 10 minutes of data has negligible consequences for the final result. We 
fit the longer-term variation with a 3rd-order polynomial and then divide it out to flatten 
the time series. We find this to be the minimum order necessary to reproduce the observed 
drifts while leaving the periodic variations largely unaffected. 

Observations with the InfraRed Array Camera (IRAC) indicate that the Spit zer Space 



Teles cope's pointing oscillation behaves rather more like a sawtooth than a sinusoid ([Morales- Calderon et al 



20061 ). We indeed find that an asymmetric sawtooth best fits the cyclic time behavior in the 
data, and that the shape of the sawtooth is fairly consistent from one slit position to the 
next. We find that the phase and asymmetry of the flux modulation (which depend on both 
the telescope motion and the distribution of imperfectly calibrated pixels) are essentially 
constant for all wavelength bins (modulo a 180° phase flip), and we adopt a fixed modula- 
tion pattern, with only the amplitude of the oscillation being a free parameter. The sawtooth 
that best fits the data for dither position 1 is shown over plotted on the drift-corrected data 
in the middle panel of Figure 1. 

After dividing the data by the best fitting periodic functi o n, we fit the data using the 



ephemerides and the secondary eclipse profile of iBakos et al. I (l2006bl ). The uncertainty in 
the predicted time of the center of the secondary eclipse is 5.4 minutes assuming a circular 
orbit. T he ephemerides we h ave adopted are consistent with updated values recently pub- 



lished by lWinn et al. I (120071 ). The timing uncertainty has been reduced to 0.5 minutes and 
is not a significant source of uncertainty in our analysis. In fitting the light curve, the two 
free parameters are the mean flux level outside of eclipse and the scaling of Bakos et al.'s 
light curve required to match the in-eclipse data. This allows us to use data taken during 
ingress and egress as well as during totality. We optimize the fit by varying the amplitude 
of the periodic function and minimizing x 2 during the in-eclipse portion of HD 189733b's 
orbit. Due to its greater extent, the mean level of the time series outside of eclipse is largely 
unaffected by changes in the periodic function, and can be independently measured to high 
precision. 

To constrain the non-periodic behavior of pixels in each wavelength bin over the course 
of six hours, we fir st divide the pre-correction, in-eclipse data with the light curve of 



Bakos et al. I (j2006bl ). We find that simply fitting to all data with t < -2500 s and t > 2500 
s is not sufficient to constrain long term drifts during eclipse and adds noticeably to the 
noise in the final spectrum. Bakos et al.'s (2006) light curve is initially scaled so that the 
upward correction at eclipse center is 0.5%. We then iterate on the final solution for each 
wavelength bin by rescaling the correction based on the previous flux estimate. The solution 
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converges rapidly in all cases, and we halt the calculation when the difference in computed 
eclipse depth between successive iterations becomes less than 0.002%. Changing the initial 
scaling of the light curve by a factor of two up or down has no effect on the final spectrum. 

Given a perfectly stable detector and platform, we would normally expect to be domi- 
nated by photon noise over the entire wavelength range for a star as bright as HD 189733. Af- 
ter making the corrections described above, the RMS spread of the individual measurements 
about the best fit light curve are indeed within 20% of the expected photon noise limit over 
most of the wavelength range. The RMS spread in the individual measurements exceeds what 
one would expect just from photon noise by up to 50% in the regions 9 — lO/mi, 11.5 — 12.5/mi, 
and 13.7 — 14.5/xm for both the Oct 21st and Nov 21st observations. Conversely, the regions 
8 — 9/im, 10 — 11/iin, and 12.5 — 13.5/im show an RMS spread that is consistent with pure 
photon statistics. This suggests that there are remaining, unmodeled temporal variations on 
the order of 0.5%, due perhaps to fringe motions and subpixel variations in responsivity. 

In Figure 2 we show the HD 189733b/HD 189733 flux ratios for each of the three 
separately analyzed data sets. The plotted uncertainties reflect only the light curve fitting 
errors and are smallest for the Nov 21st observations by virtue of their greater number and 
continuous coverage through secondary eclipse. 



4. Discussion 

HD 189733b's continuum is well measured and, multiplying the relative fluxes in Figure 
2 by the observed spectrum of HD 189733, we obtain absolute fluxes ranging from 1.8 mJy 
at 7.8/im to 0.5 mJy at 14.7/xm . Our measurement is consistent with the 660 /xJy integrated 



flux measurement at 16/mi by iDeming et al. I (120061 ). Comparing the three spectra to one 



another, the signal-to-noise ratio per resolution element ranges from about 15 at 10 microns 
to about 6 at 14.5 microns. 

The flux ratios in Figure 2 are essentially flat; the emission spectrum of HD 189733b 
over the wavelength range 7.5 to 14.7/xm is evidently consistent with that of a blackbody. 
The integrated flux ratios are 0.50%, 0.51%, and 0.46%, respectively, yielding a mean of 
0.49 ± 0.02%. Though the mean relative flux is in good agreement, the spectrum taken at 
dither position 1 appears somewhat problematic, with a number of peaks and troughs that 
are not seen at the other slit positions. Of the three, the spectrum taken at slit center is 
clearly the highest quality measurement and demonstrates the utility of keeping the telescope 
as motionless as possible for observations of this type. While there are interesting similarities 
between the spectra taken at dither position 2 and at slit center, the differences are sufficiently 
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large and numerous that we cannot identify individual absorption or emission features with 
any confidence. 

In Figure 3 we compare the relative flux spectrum me asured at slit center with a model 
of HD 189733b from lBurrows. Sudarsky. fc Hubeny I (120061 ). The model shown is for superior 
conjunction, with no clouds and 50% energy redistribution to the night side of the planet. 
The mean levels for the model and the observations are quite similar. Redward of 10/im, 
where the flux ratio is predicted to be fairly flat, the model predicts a mean relative flux 
which is about 20% higher than that of the observations. Near superior conjunction, less 
atmospheric energy redistri bution would push the re lative planet /star flux still higher. In 



contrast to the findings of lHarrington et al. I (120061 ) for v And b, it appears that strong 



day/night differences on HD 189733b are disfavored by the data. 

Blueward of 8.2 microns, the model predicts a quite significant drop in the relative flux 
due to water absorption. In contrast, the observed flux ratio s continue unabated to the blue 



end o f the spectrum. This difference is significant; scaling iBurrows. Sudarsky. fc Hubeny 



(j2006|) 's model to match the data beyond 9.5^um and summing over all bins blueward of 8.2 
microns, the model lies ~ 9a below the data. The measured flux at the blue end is not the 
product of a few bad measurements, but is reflected separately in each of the three individual 
spectra, and two independent analyses. It is also seen in simple subtractions of combined 
out-of-eclipse and in-eclipse spectra. 

A reduced planet-to-star ratio blueward of lO^tm due to water opacity is a near un i- 
versal prediction in published models of hot Jupiters to date (e.g. iFortnev et al.l (120051 ); 
Barman et al.l (120051 ) ; ISeager et al. I (120051 ) ; lBurrows. Sudarsky. fc Hubeny I (120061 )). and hence 
the absence of this signature from our data is intriguing. Simply reducing the amount of 
water in the models may not, by itself, bring them into agreement with the data. For ex- 
ample, our observations also appear to disagr ee with the low H 2 Q , high carbon to oxygen 



ratio spectrum computed for HD 209458b by ISeager et al 
drop in the flux b 



(120051 1 . 



that iFortney et al. 



In this case a similar 

ueward of 9 microns is expected due to methane absorption. We note 



(2006b) have recen tly considered dynamical models of hot Jupiter at- 



mospheres ( ICooper k. Showman 1 120061 ) to predict emergent spectra as a function of orbital 
phase, and have shown that spectral features in day-side spectra may be suppressed due to 
an isothermal pressure-temperature profile. A more thorough analysis of the consequences 
of this spectrum for models of giant planet atmospheres is forthcoming. 

The detection of HD 189733b's spectral continuum, accomplished with only 12 hours 
of Spitzer time, is a remarkable demonstration of the capability and utility of the Spitzer 
Space Telescope and the Infrared Spectrograph. Integrating over many eclipses, we can 
expect substantial increases in spectral signal-to-noise ratio, broader wavelength coverage, 
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and the detection of atmospheric constituents for both HD 189733b and other transiting 
extrasolar giant planets. While Spitzer's cryogen lasts, we can look forward to an exciting 
era of quantitative, spectral characterization of extrasolar planets, and ultimately a deeper 
understanding of planetary atmospheres under strong irradiation. 

We are grateful to an anonymous referee for several suggestions which significantly 
improved the analysis and presentation of this work. This work is based on observations 
made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology under a contract with NASA. Support for this work was 
provided by NASA through an award issued by JPL/Caltech. This study was supported in 
part by NASA grant NNG04GL22G and through the NASA Astrobiology Institute under 
Cooperative Agreement No. CAN-02-OSS-02 issued through the Office of Space Science. 
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Fig. 1. — The best-fit solution for the planetary flux in the wavelength bin centered on 7.93 
microns, showing from bottom to top the corrections applied. The flux measurements are 
shown as a function of time from secondary eclipse. The upper and lower data sets have been 
vertically offset for clarity. The filled circles in the bottom plot show the uncorrected flux 
as extracted using SPICE. The asterisks highlight data points actually used to (1) constrain 
the long term drift, (2) fit the amplitude of the periodic function, and (3) fit the eclipse 
depth. To better constrain (1) and (2), the asterisked da ta have been correcte d upwards 
during secondary eclipse using the predicted light curve of iBakos et al. I (j2006bl ). The first 
10 minutes of data are not used as they are generally affected by rapid latent charge buildup. 
The middle set of data points shows the effect of dividing by a 3rd-order polynomial, and 
over plotted is the sawtooth function that best fits the pointing-induced oscillations. The 
upper set of data points have been corrected for long term and periodic behavior, and the 
best fitting light curve is shown over plotted. The filled squares are each a mean of 30 
successive data points and serve simply as an aid in visualizing the quality of the fit. 
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Fig. 2.— The flux of HD 189733b relative to that of HD 189733, as measured at three dif- 
ferent slit positions. The spectra are unsmoothed, and the plotted uncertainties correspond 
to the formal fitting errors. 
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Fig. 3. — Compari son of the flux ratios measured a t slit c enter (filled circles) with a model of 
HD 189733b from iBurrows. Sudarsky. fc Hubeny I (120061 ). The model shown is for superior 
conjunction, with no clouds and 50% energy redist ribution to the night side of the planet . 
The upper panel shows the model as published by IBurrows. Sudarsky. fc Hubeny I (120061 ). 
while in the lower panel the model has been scaled to match the data at A > 9.5/zm. 
Intriguingly, the observed spectrum does not show the expected decrease in relative flux at 
short wavelengths due to increasing water opacity. 



